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(Abstract) The theoretical sufficient conditions for the permanence of discrete time models for competitive ecosystems with two species in patchy 
environment without or with dispersal corridor are obtained. Based on these sufficient conditions, the effect of dispersal corridor on the permanence 
of the competitive ecosystems is analyzed. Numerical simulation is given, which confirms the theoretical results. 
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1. INTRODUCTION 

Due to spatial heterogeneity and the spread of human activities, 
some habitats of biological species have been separated in 
isolated patches [1]. Dispersal corridors are established among 
these patches in order to meet the needs of certain species for 
resources, escaping for the predators or other competitors, 
exchanging the gene flow and other aspects to maintain the 
biodiversity [2-4]. At least, one purpose of the dispersal corridors 
is to avoid the extinction of the biological species in these 
ecosystems [5]. There are studies on the permanence or extinction 
of the biological species in patchy environment via theoretical 
analysis of mathematical models [1-3,5-7]. In these theoretical 
studies, many permanent or extinct conditions for the models are 
obtained, however, the effects of the dispersal corridors on the 
permanence of the ecosystems have not been clearly revealed 
through these theoretical analysis. 

Permanence is one of the basic concepts for bio-mathematical 
models on ecosystems, which is corresponding to the coexistence 
of biological species in ecosystems in large time scales [8]. In this 
paper, we first establish a discrete time model with two 
competitive species in patchy environment without dispersal 
corridor, and obtain sufficient conditions for the permanence of 
the model. Followed, a discrete time model with the same two 
competitive species in patchy environment but with a dispersal 
corridor is established, also sufficient conditions for the 
permanence of this model are obtained. After these, the sufficient 
conditions are compared and analyzed, the effects of the dispersal 
corridor are then revealed. That is, the effects of the dispersal 
corridors on the permanence of the competitive ecosystems are 
explored through theoretical analysis for the mathematical 
models. 

It is well known that numerical simulation can be carried 
out to verify the theoretical results are whether well-posed or 
not. Two numerical examples are given in this paper. One 
numerical example shows that one of the biological species in 
one patch will extinct when the parameters of the model 



without dispersal corridor are specified as some certain values 
while the theoretical permanent conditions are not satisfied, 
and the other numerical example shows that both the two 
species are permanent when dispersal corridor is introduced 
while other parameters are the same as the first example and 
the theoretical permanent conditions with dispersal corridor 
are satisfied. 

The paper is organized as follows: the theoretical 
permanent conditions for the competitive ecosystems with 
two species in patchy environment without or with dispersal 
corridor are obtained in section 2, respectively. In section 3, 
these sufficient conditions are compared and analyzed, the 
effects of the dispersal corridor are revealed. Numerical 
simulation is carried out in section 4. 

2. PERMANENT RESULTS 

2.1. Permanence without Dispersal Corridor 

Suppose that there are two competitive biological species 
living in the same habitats, denote as x and y , respectively. 
Their habitats are separated in two patches, patch 1 and 2. 
Denote x t (ri) (correspondingly, y,(»)) as the density of the 
species x ( y ) at the n' h time step. It is supposed that the 
increasing of the densities of the species x ( y ) in these two 
patches is Logistic [9], we introduce the following model 
x,(n + l) = x,(w)exp(r u -a u ^(n) -fc u y, («)), 

x 2 (n + 1) = x 2 («)exp(r 12 — a n x 2 (n) - b 12 y 2 (n)), , . . 

(1) 

y 1 (n + l) = y 1 (n)exp(r 21 -a 21 yj(w)-fc 21 ^(w)), 
y 2 (n + l) = y 2 («)exp(r 22 — a 22 y 2 (n) — b 22 x 2 (n)). 
In model (1), ( i = 1,2, j = 1, 2 ) is the intrinsic increasing 
rate of species i in patch j; a xj ( j = 1,2 ) is the intra- specific 
competitive strength of species x in patch j ; a z . ( j = 1, 2 ) is 
the intra-specific competitive strength of species y in patch j ; 
t>ij ( 7 = 1, 2 ) is the inter-specific competitive strength of 



IJICS Volume 1 , Issue 2 May 201 2 PP. 11-15 www.iji-cs.org © Science and Engineering Publishing Company 

- 1 1 - 



International Journal of Information and Computer Science 



IJICS 



species y with x in patch j ; b 2j ( j = 1, 2 ) is the 
inter-specific competitive strength of species x with y in 
patch j . All the parameters in model ( 1) are positive constants. 
The initial values of model (1) are 

x, (0) > 0, x 2 (0) > 0, y 1 (P)> 0, y 2 (0) > • (2) 
It is obvious that the solutions of model (1) with initial 
values (2) are positive that is consistent to the biological 
background. 

Next we discuss the permanence of model (1) with initial 
values (2). we give the following lemma which can be easily 
obtained from Lemmas 1 and 2 in [10]. 

Lemma 2.1. Suppose that r, a are both positive constants, 

{x(n)}™ =1 is a positive sequence. If 

x(n + 1) < x(n) exp(r - ax(n)) , 

then 

limsupx(n) < — exp(r-l); 

n— >co Q. 

If {x(n)}™ =1 satisfies 

x(n + 1) > x(n) exp(r - ax(n)) 
and limsupx(n) < k ( k is a positive constant), then 



limsupy 2 (w) < A 22 . 



(11) 



provided that 



liminf x(n) >— exp(r-ak) 



-£>1. 



For the simplicity, we denote 

Ai = — ex P(li "I)' A 2 = — exp(r 12 -1), 



1 



*12 
1 



(3) 



Ai = — exp(r 21 -1), A 22 = — exp(r 22 -1). 
a,. a. 



*21 "22 

The next theorem gives the sufficient conditions for the 
permanence of model (1) with initial values (2). 
Theorem 2.2. If 

r u -V*2i>0, (4) 
r 12 -fc 12 A, 2 >0, (5) 
r 21 -b 21 A n >0, (6) 
r 22 -fo 22 A 12 >0, (7) 

then model (1) is permanent with initial values (2). 

Proof. From the first equation of model (1), 

Xj(n + 1) < Xj(tt)exp[r u -a^iri)] . 

Hence, from Lemma 2.1, 



Similarly, 



limsupxj(n) < — exp(r n -1) = A,, . 



limsupx 2 (n) < A 12 , 

rt->=o 

limsupy,(n) < Aj, , 



(8) 

(9) 
(10) 



From (10), there exists a positive integer ^ such that for any 
s > one has for " > ^ 

y 1 («)<A 21 +;r. (12) 
From (12) and the second equation of model (1), 

Xj(n + 1) > x 1 (n)exp[r 11 -b^A,, +s)-a n x i (n)] , 
further, from Lemma 2.1, 

liminf x { (n) >— — ^"^Ai +g ) eX p[ r]i -b^A^ +e)-a n A u ) , 



provided that 



Denote 



r n -b n (A 2i +s) 



A U >1. 



(13) 
(14) 



r n Ail 



exp(r 11 -fo 11 A 21 -a 11 A 11 ), 



> is consequent from (4) and note that £ is arbitrarily 
given, from (13), (14), 

liminf x^ri) > >0 , 

»-»co 

if 

fl, 



holds. Again from (4), we get 



" Ai>i, 



(15) 
(16) 



a u A _ exp(r n -l) ^ r n > ^ 
r n -b n A 2l r n -b n A 21 r n -b n A 21 
That is, (16) is true if (4) holds. Moreover, 

< < < ^ < — exp(r u -1) = A u , 

since r n -b n A 2l -a u A n < from (16). As a whole, 

< < lim inf x, (n) < lim sup jc, (n) < A n . (17) 

Similarly, there exists positive constants m' 2 ' , m ( 2 , m 2 2) such 
that 

< m, <2) < liminf x 2 {n) < limsupx 2 (n) < A 12 , (18) 
0<m™ < liminf y,(n) < limsupy,(n) < A 21 , (19) 
< mf ' < liminf y 2 (n) < lim sup y 2 (n) < A 22 . (20) 
Thus completes the proof of Theorem 2.2. 

2.2 Permanence with Dispersal Corridor 

In this part, we suppose that the two isolated patches in section 
2.1 is now introduced one dispersal corridor. The species x 
can migrate between patch 1 and 2 with this corridor freely, 
but the species y can not disperse between these two patches, 
in this case, model (1) is rewritten as 
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+ = x,(n)exp(r H - a n x { (n) -b^^n) - d l2 x { {n) + d 2l x 2 (n)), 
x 2 (n + l) = x 2 (n)exp(r u -a l2 x 2 (n)-b l2 y 2 (n)-d 2l x 2 (n) + d l2 x l (n)), 
y, (n + 1) = y l («) exp(r 21 - a 21 y l (n) - b 2X x x (n) - D 12 x, («) + D 2 ,x, (n)), 
y 2 (n + 1) = y 2 (n) exp(r 22 - a 22 y 2 (n) - £ 22 x 2 (n)-D 2l x 2 («) + D 12 x, (n)). 

(21) 

In model (21), the positive constant aL is the intra- specific 
competition strength due to the migration of species x from 
patch i to j ( ;', j = 1, 2 ), and the positive constant Z>. is the 

inter-specific competition strength due to the migration of 
species x from patch ; to j ( i, j = 1,2 ). It is obvious that 

model (21) is the same as model (1) if d v = Z> = , ;', j = 1, 2 . 



exists a subsequence {x 2 (n j )}'° =1 of {x 2 (n)} , that satisfies 



Theorem 2.3. If 



(22) 
(23) 
(24) 
(25) 

then model (21) is permanent with initial values (2), where 
exp(r u 



I1-M21 >0 > 
r l2 -b l2 A 22 >0, 

r 21 -(fc 21 +D 12 )B n >0, 

r 22 _ 

(£ 22 +Z) 21 )fi 12 >0, 



!) d exp(r 12 -l) 

",«12=- 



fl n +J 12 ' " a 12 +rf 21 
Proof. Considering the auxiliary system 

w(n + l) = M(«)exp[r u -(a u + d l2 )u(n) - b u v(n)], 
[v(n + 1) = v(n) exp[r 21 - (b 2] + D l2 )u(n) - a 21 v(n)], 
with initial values w(0) > 0, v(0) > . From the proof of 
Theorem 2.2, this system is permanent if (22) and (24) hold. 
Hence, there exists positive constants m, (1) and ni 2 such that 

(i) 



liminf u(n) > m? 



liminf v(n) > m 2 

n— >co n— >co 

But for any positive integer n , 

x, (n + 1) > u(n + 1), V[ (n + 1) > v(n + 1), 

hence, 

liminf x,(n) > mj 11 , liminf Vj(m) > m ( 2 . 

H->co n->co 

Similarly, 

liminf x 2 (n) > mj 2 ', liminf y 2 (ri) > m 



(2) 



(26) 
(27) 



where mf 11 and m 2 ' are positive constants. 

Next we prove the following to obtain the permanence of 
model (21) with initial values (2): 

lim sup Xj («) < Mj (1) , lim sup x 2 («) < Mj <2) , 

(28) 

limsup>'! (n) <M™,limsupy 2 (n) <M { 2 \ 

where M, 0) (i,y=l,2) is a positive constant. 
Suppose the contrary, there exists no M, (1) such that 
limsupX[(«) <Mf\ 

therefore, there exists a subsequence {x^nj}"^ of {Xj(n)} , 
that satisfies 

lim x, (n t ) = +00 . 

k— >co 

Further, from the second equation and Lemma 2.1, there 



lim x 2 (n ; ) = +00 . 



(29) 



Multiplying the first and the second equation of (21), 

Xj (n + l)x 2 (n + 1) < x 1 (n)x 2 («)exp[r 11 + r 12 a 11 x 1 (n)-a 12 x 2 («)] 

Note that the function f(x, y) = xy exp(r - ax - by) is 
bounded in the first quadrant, one has 

lim sup X[ (n + l)x 2 (n + 1) 

exists finite. Note (29), limx,(n )=0, which contradicts 

/-MO 

with the first inequality of (26). This contradiction shows that 
the first inequality of (28) holds. The other three inequalities 
of (28) can be similarly proven. Therefore, model (21) with 
initial values (2) is permanent from (26), (27) and (28). The 
proof of Theorem 2.3 is completed. 



THE EFFECTS 
CORRIDOR 



OF DISPERSAL 



The sufficient conditions for the permanence of model (1) and 
(21) have been obtained in section 2, respectively. In this 
section, we compare and analyze these two groups of sufficient 
conditions, that is, conditions (4)-(7) and (22)-(25). 

The condition (4), (5) is the same as (22), (23), respectively. 
Condition (6), (7) can be written in the parameters of model (1) 
as 



and 



-exp(r u -1) > , 



r 22 - — exp(r 12 - 1) > , 



respectively. And condition (24), (25) is equivalent to 



^ + %exp(r n -l)>0, 



a n +d i: 



and 



-^W 12 

a l2 +d 2l 



-1)>0, 



(30) 



(31) 



(32) 



(33) 



respectively. Comparing (30) with (32), and (31) with (33), the 
difference is that the terms d tj and D tj (i, j =1,2 ). Recalling 
model (1) and (21), these terms reflects the intra- or inter- 
specific competition strength due to the migration of the species 
x . It can be observed that the positive term d tj in the denomi- 
nator is helpful to satisfy both the conditions, (32) and (33), 
while the positive term Z> in the numerator is a disadvantage 
both for the conditions, (32) and (33). Hence, the migration of 
the species x through the dispersal corridor is helpful for the 
survival of themselves, but the migration of the species x is a 
disadvantage for their competitors, the species y , since 
additional inter-specific competition is produced due to the 
dispersion of the species x via dispersal corridor. 
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But the 



values of 



b 2 2 +D U 



are affected 



simultaneously by the terms d tj and D„ ( i, j = 1,2 ), therefore, 

the effects of dispersal corridor on the permanence of the 
competitive ecosystems are relevant to which competition 
strength due to the migration, the intra-specific competition or 
the inter-specific competition, is stronger. 

4. NUMERICAL SIMULATION AND 
DISCUSSION 

In this section, we carry out numerical simulation to confirm 
that the theoretical results of section 2 are well-posed and 
illustrate the explanation in section 3 for the effects of dispersal 
corridor on competitive ecosystems. 

Example 4.1. Consider the following special case of model 
(1): 

x^ti + 1) = x l (n)exp(03-0.2x l (n) -O.ly^n)), 
x 2 (n + l) = x,(«)exp(0.52-0.35x,(«) -0.l5y 2 (n)), (^4) 
y l (n + 1) = y l (n) exp(0.2 - 03y t (n) - 0. 15jc, (m)), 
y 2 (n + 1) = y 2 (n) exp(0.4 - 0.25y 2 (n) - 0.2x 2 («)). 
In model (34), we have 

r n =0.1502, r n -b l2 A 22 =0.1907, 

r n -KAx = -0.1724, r 22 -b 22 A 12 =0.0464, 
that is, condition (6) is not satisfied. Hence, model (34) is not 
permanent by Theorem 2.2, and it can be observed from Figure 
1 that the species y in patch 1 will extinct. 




ICO 130 200 HO JOB 350 400 450 500 



In model (35), dispersal corridor is introduced for the species 
x , and the competition strengths produced by the migration are 
d X2 = 0.19, £>,, = 0.001 and d 2l = D 21 = 0, while the other 

parameters are the same as model (34). The conditions (22)-(25) 
are 



r n =0.1502, r n 



-b X2 A 22 =0.1907, 



r 21 - (b 2l + D n )B U = 0.0077, r 22 - (b 22 + D 2l )B l2 = 0.0464 . 

That is, the conditions (22)-(25) are all satisfied, hence, 
model (35) is permanent by Theorem 2.3, which is illustrated 
by Figure 2. 

Model (34) together with model (35) shows the sufficient 
conditions of Theorem 2.2 and 2.3 are well-posed. Moreover, 
these two examples imply that the dispersal corridor is helpful 
for the permanence of the competitive ecosystems if the 
migration produces weaker inter-specific competition and 
decreases some stronger intra-specific competition at the same 
time, that is corresponding to the explanation of the last 
paragraph in section 3. 

The effect of the dispersal corridor on competitive 
ecosystems is discussed with two species and two patches in 
this paper, It should be pointed out that Theorem 2.2 and 2.3 
can be extended to competitive ecosystems with more species 
and patches. 




Figure 1. The variations of the population densities of model (34). The vertical 
axis is the population densities and the horizontal axis is the time. The red line 
is the population density of the species x in patch 1, the black line is the 
population density of the species x in patch 2, the green line is the population 
density of the species y in patch 1 , and the blue line is the population density 
of the species y in patch 2. The species y in patch 1 will extinct, and model (34) 
is not permanent. 

Example 4.2. Consider the following special case of model 
(21): 

x, (n + 1) = x, (n) exp(0.3 - 0.2x, («) - 0. 1 v, (n) - 0. 19x, («)), 
x 2 (n + 1) = x 2 (n) exp(0.52 - 0.35x 2 (n) - 0. 15 y 2 (n)), ( 35 ) 
y, (n + 1) = y, (n) exp(0.2 - 0.3y, (n) - 0. 1 5x, («) - O.OOlx, («)), 
y 2 (n + 1) = y, (n) exp(0.4 - 0.25 y 2 (n) - 0.2x 2 (n)). 
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Figure 2. The variations of the population densities of model (35) that is 
permanent. The other aspects are the same as Figure 1 . 
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